The role of disorder symmetry in a dirty superconductor: a Bogoliubov-deGennes 

(BdG) study 
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We study the effect of disorder symmetry in a disordered s-wave superconductor. We begin 
with an attractive Hubbard model on a two-dimensional square lattice in the presence of diagonal- 
and off-diagonal disorder. The model is studied numerically using Bogoliubov-deGennes approach, 
which proceeds via the decoupling of attractive Hubbard term in the pairing channel. We find that 
the off-diagonal disorder is much more effecient in destroying superconducting order. A detailed 
analysis brings out very distinct qualitative pictures in the two cases: diagonal disorder leads to 
formation of small scale isolated superconducting islands separated by large region of normal metal, 
whereas off-diagonal disorder leads to a connected network of superconducting region. Naively, it 
seems that the connected network of superconducting regions should have a large average value for 
superconducting order parameter, but this does not hold true for large disorder, as we demonstrate 
explicitly. Our qualititive picture is supported by real-space data on local order parameters. 



Introduction: The behavior of a fermion system in the 
presence of disorder and strong inter-particle interactions 
has been a challenging question in condensed matter 
physics for many decadesii"— Anderson's seminal discov- 
ery that the wave nature of electrons in a quantum me- 
chanical context can cause the electrons to completely 
localize in a disordered system gave rise to the field of 
Anderson localization.i^ On the other hand, in the clean 
system, the Coulomb repulsion between the electrons 
themselves can also cause the electrons to localize, a phe- 
nomenon known as Mott localization,— the most famous 
example of such materials being the undoped parent com- 
pounds of the high-Tc Cuprate superconductors<i 

A well-studied problem in condensed matter physics 
has been the properties of a dirty superconductor, with 
early theories by Anderson^ and Abrikosov and Gorkox*^ 
in the regime when the disorder was weak. With the 
advent of advanced experimental techniques of grow- 
ing highly inhomogenoeous films, research on the prop- 
erties of a dirty superconductor in the strong disorder 
regime has become possible.— This is also the regime 
where electron interactions start to play an important 
role, which then becomes a very challenging theoretical 
problem)^ Also, controlled studies of the related quantum 
phase transition between the superconductor and the An- 
derson Insulator (SIT) have been widely carried out. In 
the SIT, the early theoretical prediction of a universal 
conductivity (in units of e^/S)^ demarcating the criti- 
cal point between the two phases generated a plethora of 
theoretical and experimental work in this fieldf^^rJ^ Re- 
cently, with the observation of Anderson localization in 
disordered optical lattices and the possibility of tuning 
the sign of the interaction between the atomic species, 
a new avenue has opened up to investigate the interplay 
of disorder and interaction among both fermionic and 
bosonic atoms J^"— 

In the presence of attractive interactions between 



fermions, disordered systems can become superconduct- 
ing (for charged electrons) or superfluid (for neutral 
atoms). Below a critical disorder strength in d = 2, the 
two-particle states may sustain dissipationless fiow even 
though all single-particle states are localized due to An- 
derson localization. The combination of attractive inter- 
actions and disorder can lead to a rich variety of phases 
and new physics, including the BCS-BEC crossover, var- 
ious kinds of glassy phases and superfluids (superconduc- 
tors). 

A related very important aspect of physical systems 
with disorder is the symmetry of the Hamiltonian gov- 
erning them. A complete classification of the symme- 
tries of the single-particle Hamiltonians governing a dis- 
ordered system was recently achieved.—"— Depending on 
the presence (or absence) of time-reversal symmetry (T) , 
spin-rotation {SU{2)) symmetry and particle-hole sym- 
metry (S) in the Hamiltonian, the physics of the disor- 
dered system can be markedly different, particularly with 
respect to Anderson localization. For example, it is now 
known that weak localization does not occur in the chiral 
symmetry class in the conventional sense, but follows a 
topological route which then gives way to strong local- 
ization for stronger disorder)^ 

In traditional solid state systems, intrinsic disorder 
typically originates from random locations of dopant 
ions. This invariably affects both the on-site potential 
and the hopping parameters in the corresponding model 
Hamiltonian. Depending on the details of the structure 
and ions involved, this may affect one or the other more 
strongly and therefore can be modelled via either a ran- 
domness in on-site potentials (diagonal disorder) or hop- 
ping strengths (off-diagonal disorder) . In addition, disor- 
dered optical lattices provide another promising direction 
to investigate the infiuence of symmetries on the local- 
ization problem. 

The SIT has been investigated in considerable detail 



2 



in the literature with orthogonal disorder i.e., on-site 
or diagonal disorder with and without the spin-rotation 
(SU(2)) symmetry.—"— In this paper, we investigate the 
influence of off-diagonal disorder in the physics of a dirty 
superconductor, which belongs, in the classification of Al- 
tland and Zirnbauer— in the chiral-BdG symmetry class, 
and compare it to similar physics in the usually studied 
BdG-orthogonal symmetry class. Even though we do not 
explicitly investigate the physics of the critical point it- 
self, the effect of imminent localization is always present 
in the physics of a dirty superconductor, and the nature 
and the properties of the inhomogeneities in the super- 
conducting state itself is a very important indicator of 
the subsequent physics of localization. Elucidating how 
the superconducting state responds to the nature of the 
disorder potential will be the main goal of this paper. 

Model and Method: The starting point of our investi- 
gation is the attractive Anderson-Hubbard model on a 
square lattice, given by the Hamiltonian 



sp 



(1) 



where Hgp is the single-particle part and Hint is the in- 
teraction part. Hsp is given by 



(2) 



Here Cia (c^o-^ annihilates (creates) an electron at site Ri 

with spin-projection c, {ij) implies that Ri and Rj are 
nearest neighbors, /i is the chemical potential and tij is 
the hopping matrix clement between sites Ri and Rj . We 
write tij (= tji) as 



U 



t + 6, 



(3) 



where Sij is drawn randomly from a box probability 
distribution P (Sij) = -^0 - \Sij\), & (x) being the 
Heavisidc 0-function. Similarly, the on-site diagonal 
disorder term is denoted by the random site-energy e^, 
which is drawn from a box-distribution given by P (si) = 
yQ (y — l^il)- The attractive Hubbard interaction term 
Hint is given by 



PIE 



ni-fUii, 



(4) 



where 



is the number operator at site Ri 



and spin-projection cr. We set t = 1 as the basic scale 
for energy, therefore all energy parameters are in units 
of t. This leaves four independent energy-scales in the 
model: the interaction strength \U\, the off-diagonal dis- 
order bandwidth Vt , the diagonal disorder bandwidth V 
and the temperature T. If Vt = 0, the model reduces to 
the class of disorder Hamiltonians with orthogonal sym- 
metry. If y = 0, it reduces to the class of disorder Hamil- 
tonians with chiral symmetry. 

At half-filling, the clean (i.e., no disorder) attractive 
Hubbard model (in d = 2) has a phase diagram in which 



the charge density wave (CDW) state and the supercon- 
ducting state (SC) are exactly degenerate in energy and 
co-exist at T = 0. As one dopes the system away from 
half-filling, the CDW state gradually disappears, and the 
system makes a Berezinskii-Kosterlitz-Thouless (BKT) 
transition to a conventional s-wave SC ground state.— 

A discussion of the effect of non-magnetic disorder on 
superconductivity (in particular, Tc) begins with Ander- 
son's theorem^ which states that, at least for weak impu- 
rity scattering, Tc is not affected by disorder. In the con- 
text of mean-field BdG formalism with diagonal disorder, 
Ghosal et found that it is impossible to destroy 

superconductivity through amplitude fluctuations by in- 
creasing disorder (an extreme case of Anderson's theo- 
rem), and in fact, there is a disorder strength for which 
the spectral gap (the superconducting gap) shows a non- 
zero minimum, beyond which the gap actually increases 
with disorder strength, even though the spatially aver- 
aged SC order parameter monotonically vanishes. Con- 
ventional wisdom suggests that the impurity states will 
gradually fill in the superconducting gap and destroy SC. 
However, as we shall see, this intuitive picture is not al- 
ways true, and in fact, depends on the type of disorder. 

We revisit this scenario where the disorder resides on 
the bonds (ofT-diagonal/bond/hopping disorder) rather 
than on the sites (diagonal/sitc/potential disorder). The 
BdG equations are derived as usual, by decoupling the 
attractive Hubbard interaction in the pairing channel. 
This leads to the effective BdG Hamiltonian)^ 



{ij),(y if ia 



(5) 



where liia = \i^\U\ (ni._o-). i^off is diagonalized via the 
Bogoliubov transformation: 



n 



t * 
„t „* 



(6) 



where 7„cr (tIct) ^'^^ annihilation(creation) operators for 
the Bogoliubov quasiparticles in the single-particle state 
(n, cr) and satisfy canonical anticommutation relations for 
fermions. The local superconducting order parameter 
and particle densities arc obtained self-consistently via 

A, = -|C/|(c,tc.^) 

= -\U\Y,!{E,:)u.^vl^ 



71 } "-n 



n 



|2 



(7) 



where is measured from the global chemical poten- 
tial, fj, and f {En) is the corresponding Fermi function. 
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The original Hamiltonian in Eqn.[T]and the effective BdG 
Hamiltonian in Eqn. [S] are manifestly invariant under T 
and SU{2). li V = (but Vt ^ 0), the energy eigenval- 
ues also appear in ± pairs about the chemical potential, 
thus exhibiting chiral symmetry. In the superconducting 
state, the disordered BdG Hamiltonian thus belongs to 
the BdG chiral class. 

All calculations are performed in the grand canoni- 
cal ensemble scheme by keeping the chemical potential 
H constant and allowing the average electronic density n 
to vary freely. This approach is conceptually accurate 
since no adjustments are required in order to achieve 
a target average density at the end of the calculation. 
The disadvantage is that we need to perform calculations 
for a range of /i values in order to achieve the desired 
n. However, this approach naturally provides us results 
for density dependence as well. The consistency loop is 
started by picking a random set for the local order pa- 
rameters Ai. Subsequently, numerical diagonalization is 
carried out and the eigenvectors and eigenvalues are used 
to compute the new set of A^, and so on until the order 
parameters and local densities converge to a desired ac- 
curacy. In this work we use a convergence accuracy of 
10""' for both the and the charge density, i.e., we stop 
the self-consistency loop when the value in (fc + 1)*^ step 
differs from the one in k^^ step by less than 10^^. 

Results: 

We begin our discussion by showing the behavior of the 
spatially averaged superconducting order parameter Aop 
with disorder strength, V, in Fig. [ija) and (b). Here, 
Aop is defined as 



A, 



(8) 



where Ns denotes the number of lattice sites and the 
overbar denotes an average over different realizations of 
quenched disorder. 

In Fig. [TJ panel (a) shows the behavior of a system 
with diagonal (site-) disorder and panel (b) shows the 
behaviour of a system with off-diagonal (bond-) disorder. 
In both cases, we present results for systems where the 
linear lattice size takes values ranging from L = 24 to 
L — 64. It is clear from Fig. [1] that Aop has very little 
system-size dependence and has, in fact, converged very 
well already for the system sizes L = 24 and L = 32. 
We have ascertained that, in this and all our subsequent 
results, all quantities have converged by increasing the 
system size. 

In Fig. UKa) and (b), there is already a clear qualita- 
tive difference in the behaviour of the superconducting 
order parameter for the two kinds of disorder as disorder 
strength is gradually increased. For the diagonal disorder 
case (Fig.[TJa)), the order parameter is seen to gradually 
converge to a non-zero value as disorder strength V/t is 
increased. Indeed, even when the disorder strength is as 
large as the kinetic bandwidth Wkin = 8i, significant non- 
zero superconducting order parameter is found, when av- 
eraged over the entire system. This was understood in 




FIG. 1: (Color online) Dependence of the average supercon- 
ducting order parameter, (a)-(b), and the spectral gap, (c)- 
(d), on disorder strength for (a),(c) diagonal(site-) disorder 
and (b),(d) ofr-diagonal(bond-) disorder for various lattice 
sizes. The plots are shown for \U\ /t — 3.0 and n ~ 0.88. The 
disorder realization average varies from 20 copies for L = 24 
to 4 copies for L — 64. The inset in panel (b) compares 
the values of Aop between diagonal and off-diagonal disor- 
der at large disorder bandwidths. Note that the value for 
off-diagonal disorder is smaller by an order of magnitude. 



terms of a phase separation of the system between super- 
conducting islands (with large values of A^) separated by 
insulating normal regions with very small A^i^ 

On the other hand, the figure is very different for off- 
diagonal disorder. As disorder strength is increased in 
Fig. [Ijb), the order parameter supresses rapidly. In or- 
der to clarify the difference, we plot the comparison on a 
logarithmic scale in an inset in Fig.[T](b). It is clear from 
the inset that at large disorder strengths Aop is almost 
an order of magnitude smaller for the bond-disorder. In 
other words, at the level of mean superconducting or- 
der parameter, bond-disorder destroys superconductivity 
very quickly. 

The difference in behaviour of the dirty superconduc- 
tor for the two kinds of disorder becomes even more clear 
if we consider how the superconducting spectral gap re- 
sponds to a change in the disorder strength. In Fig. [TJc) 
and (d), we plot the spectral gap -Egap defined as 



E — E — E 



(9) 



Here E'oxc is the energy of the first excited state and Eg 
is the ground state energy. The overbar, once again, de- 
notes averaging over various realizations of quenched dis- 
order. Our results for the diagonal disorder (panel (c) in 
Fig.[T]) compare very well with those presented in Ghosal 
et al^ This makes us confident that the behaviour of 
the superconducting gap is a physical property of the 
system, not a numerical artefact. Similar to Fig. [Ija) 
and (b), we plot the spectral gap for increasing lattice 
sizes. For off-diagonal disorder, Fig. [Ud), the spectral 
gap converges very well for all lattice sizes. Some minor 
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FIG. 2: (Color online) Disorder and density dependence of 
the average superconducting order parameter and the spec- 
tral gap for diagonal (filled symbols) and off-diagonal (open 
symbols) disorders. Results are shown for two values of U/t 
as indicated in the figure. Note that the superconducting or- 
der parameter as well as the gap show a significant density 
dependence in case of off-diagonal disorder. 



size-dpendence is seen for smaller lattices in the case of 
diagonal disorder, Fig. [TJc), even though the results are 
manifestly independent of the lattice size. 

For diagonal disorder, we confirm a counter-intuitive 
behavior first observed by Ghosal et al^^^: the spectral 
gap initially decreases with an increase in V, but after a 
minimum around ~ 4i, it monotonically increases with 
V. Ghosal et al. ascribed it to the fact that at very high 
diagonal disorder, the non-interacting eigenstates are es- 
sentially spatially local, and only those local eigenstates 
near the chemical potential participate in pairing. The 
gap equation in this local approximation implies that the 
gap depends inversely on the square of the localization 
length (Cioc), which decreases with increasing disorder. 
This causes the increase in the spectral gap with increas- 
ing disorder. 

However, for off-diagonal disorder. Fig. [IJd), the be- 
haviour of the spectral gap is strikingly different. With 
increasing Vt, i?gap monotonically decreases. For the 
largest system sizes with converged values for -Egap, it 
is seen to converge to a very small value - again reaffirm- 
ing the fact that off-diagonal disorder is more effective 
compared to diagonal disorder in destroying supercon- 
ductivity. 

In Fig. [2]pancls (a)-(d), we introduce the dependence 
of the superconducting state properties, Aop and Sgap 
on the average electron density, n. In Fig. [D^a), Aop 
is plotted against the electron density for U = 2t and 
for various values of V/t (full symbols with full lines) 
and Vt/t (open symbols with dashed lines). We use the 
same set of values for both V and Vt (from V = 0.2 



V(V,)/t V(V,)/t 

FIG. 3: (Color online) Disorder dependence of the supercon- 
ducting order parameter for different values of n, and for the 
diagonal and off-diagonal disorder. Note that the relative ef- 
fect of disorder depends on the average electronic density, n. 



to V ~ 4.0), so that comparisons become meaningful. 
For U/t = 2, Fig. [D^a), the magnitude of Aop decreases 
monotonically with V/t. At weak disorder, Aop is weakly 
density-dependent for both types of disorders. As the 
disorder bandwidth increases, the Aop for the diagonal 
disorder becomes essentially independent of n, whereas, 
the Aop for off-diagonal disorder continues to have a sig- 
nificant density dependence. Interestingly, there seems to 
be a crossover density Ucr above which bond-disorder has 
stronger superconductivity than site-disorder (in terms of 
having a larger Aop). For example, at half- filling, Aop is 
consistently larger for Vt than for the same V, for all 
disorder strengths investigated. Fig. EJb) shows that the 
same conclusion is true for U/t = 3.0, except that ricr is 
pushed towards lower densities. However, Fig. [ija) and 
(b) show that for even larger disorder strengths (for ex- 
ample, at y = Vt = 6i, n ~ 0.88 and U/t = 3.0), this 
trend will again be reversed and the order parameter, 
although small, will be larger for site-disorder than for 
bond-disorder. 

In Fig. [2] panels (c) and (d), we plot the supercon- 
ducting energy gap for two values of U, U = 2.0t and 
U = 3.0i. The non-monotonicity in the energy gap for 
site-disorder is seen to occur for all densities and for both 
the values of the interaction strength, indicating that this 
is a ubiquitous phenomenon for a dirty superconductor 
with on-site disorder. On the other hand, it is absent for 
off-diagonal disorder for any parameter regime. The dif- 
ferences in the density dependence for the two disorder 
types also show up in the superconducting energy gap. 
The plots for diagonal disorder are essentially flat for 
V = 2t (filled down triangles) for both values oiU/t. The 
corresponding plots for the off-diagonal disorder show a 
relatively strong density dependence. 
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(a) V=0.5t (b) V=t 




(c) V=2t (d) V=4t 




FIG. 4: (Color online) Real space picture of the magnitude 
of Ai for different strengths of the diagonal disorder. 



In order to emphasize the qualitative differenees be- 
tween the density dependence for the two disorder types, 
we plot Aop for four representative values of the average 
electron density in Fig.[3l panels (a) - (d). Close to half- 
filling (n = 1), the bond-disorder has larger average order 
parameter, while the site-disorder dominates for all dis- 
order strengths as we move away from half-filling. The 
difference in the density dependence points to a quali- 
tative difference between the manner in which the two 
disorder types affect the superconducting state. 

In order to investigate further the differences between 
the effects of the two disorder types, we present the real- 
space data in Fig. |4] and Fig. [5l We plot the real space 
distribution of the local order parameter for both diago- 
nal and off-diagonal disorder (Fig. |4] and Fig. [5|) and at 
the interaction strength U = 3.0t and the average elec- 
tron density n 0.88. Looking at the real space distri- 
bution of local quantities provide important clues about 
the underlying physics. In Fig. |4l for diagonal disorder, 
we clearly see the formation of disconnected supercon- 
ducting islands which has large values of the supercon- 
ducting order parameter Aop, separated by large normal 
regions with no superconducting order. As the disorder 
strength is increased, the separation of the homogeneous 
superconducting system into superconducting and nor- 
mal regions become more and more pronounced. For 
very large diagonal disorder, it is expected then that the 
system will break up into nearly microscopic regions of 
superconducting islands, interspersed by normal regions 
in between. 

In Fig. [5l for off-diagonal disorder, we see a different 
behavior of the real-space distribution of Aop (note that 
we show the same disorder bandwidths for both diagonal 
and off-diagonal disorder, so that comparisons are mean- 
ingful). For off-diagonal disorder, even for high disorder 
strengths, the distribution of the superconducting order 
parameter is much more homogeneous, not only in space 
but also in magnitude. There are very few lattice sites 
(or regions) with very high values of A^. Thus, for the 



off-diagonal disorder, superconductivity is spread more 
homogeneously but weakly over the entire system as op- 
posed to diagonal disorder. We see a similar behavior for 
all the individual disorder configurations we have inves- 
tigated. 



(a) Vj=0.5t (b) V,=t 




(c) V,=2t (d) V,=4t 




FIG. 5: (Color online) Real space picture of the magnitude 
of Ai for different strengths of the off-diagonal disorder. 

In order to make the comparison more quantitative, we 
plot the distribution, N{A), for the local order parame- 
ters in Fig. [6] N{A) measures the number density of sites 
that have the value of superconducting order parameter 
belonging in the range (A — 6/2, A -f (5/2), with S — 0.01. 
Fig. EJa) shows the distribution for diagonal disorder for 
three values of the disorder bandwidth. As the disorder 
bandwidth is increasing, there are always sites having a 
large value for the superconducting order parameter, con- 
tributing to the right tail of the distribution. In fact, the 
cut-off value for the distribution increases, as indicated 
by the horizontal arrow in the figure. On the other hand, 
the distribution for off-diagonal disorder shows that the 
superconductivity is suppressed everywhere with increas- 
ing disorder bandwidth and no sites remain with high 
values of the superconducting order parameter. Indeed, 
the cut-off value decreases rapidly in this case. This be- 
havior is consistent with the real-space pictures discussed 
above. 

Discussion: In the absence of a full analytical solution 
to the problem of a superconductor in the presence of 
strong disorder, any explanation of the different behav- 
ior of the dirty superconductor with diagonal and off- 
diagonal disorder can be at best heuristic. Typically, 
in a disordered system with on-site disorder, sites with 
the similar disorder potential reside far away from each 
other. In other words, the local potential on two neigh- 
boring sites are typically very different. This impedes the 
hopping of the electron over any macroscopic distance, 
leading to local pair formation. As disorder bandwidth 
increases, the system thus finds it advantageous to sep- 
arate into smaller and smaller microphases of supercon- 
ductor and normal regions (whether the normal regions 
are metallic or insulating will be discussed elsewhere). 
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FIG. 6: (Color online) The distribution of local amplitudes for 
the superconducting order parameters (Ai) for (a) diagonal 
and (b) off-diagonal disorder. The horizontal arrows indicate 
the direction in which the maximum value of Ai moves with 
increasing disorder bandwidth. 

However, the same is not true for ofF-diagonal disorder. 
A site which has a bond that impedes hopping because 
of very small hopping amplitude, will typically also have 
another bond which leads to a significant gain in kinetic 
energy, i.e., promotes hopping. Therefore, a connected 
network-like region appears with an almost homogeneous 
order parameter (see Fig. EJd)). Thus the pairing am- 
plitude is uniformly spread over the system even when 
the randomness in the hopping matrix element is signif- 
icant. To summarize the differences, for onsite disorder 
the order parameter gets contributions from a few islands 
with large magnitudes of the order parameter, while for 
bond disorder, the order parameter gets contributions 
from all over the system, though, with increasing disor- 
der strength, the contribution becomes uniformly weak. 
This qualitative picture is very well supported by the dis- 
tributions of the local order parameter values shown in 
Fig.m 

In fact, the qualitative difference in the effective pic- 
tures also explains the difference in the density depen- 
dence that we reported in Fig. [5] and Fig. [3] In case 
of diagonal disorder, since the microscopic regions with 
large value for the average superconducting order param- 



eter are isolated, the extra electrons find a different island 
to occupy and the system remains insensitive to the elec- 
tronic density. On the other hand, bond disorder essen- 
tially leads to a connected sub-system which can be con- 
sidered as a sub-system of a new shape and size with weak 
disorder and therefore has a density dependence simi- 
lar to that in a weakly disordered system. We propose 
qualitatively distinct pictures of how disorder destoys su- 
perconductivity in a site- vs. a bond-disordered system. 
We believe that this also has implications for how super- 
conductivity will be destroyed at the critical point be- 
tween the superconductor and the Anderson insulator. In 
agreement with previous theoretical work)^i^i2£, for the 
diagonal disorder, the non- vanishing spectral gap implies 
that even with increasing disorder, low-energy fermionic 
quasiparticles are not generated. Thus the superconduc- 
tivity will be destroyed by increasing phase fluctuations 
among the superconducting islands. The critical theory 
will be bosonic, and in c? = 2, will consist of a diffusing 
metal of Cooper pairs at the quantum critical point itself. 
On the other hand, the superconductor in the presence of 
strong off-diagonal disorder will be destroyed in a differ- 
ent fashion, with the order parameter uniformly vanish- 
ing across the system, and with the presence of fermionic 
quasiparticles at the transition (for a recent experimental 
work on the possibility of the existence of a bosonic and a 
fermionic SIT, see HoUen et al)S^ Unfortunately, in con- 
ventional solid state systems, the tuning of the type of 
disorder is experimentally very difficult. However, with 
the advent of the cold atom systems, we hope different 
scenarios with different types of disorder can be tested 
experimentally. 
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